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ABSTRACT 

We calculate numerically the radiation spectrum from relativistic electrons 
moving in small scale turbulent magnetic fields expected in high energy astro- 
physical sources. Such radiation spectrum is characterized by the strength pa- 
rameter a = Abc I S I /mc^, where Ab is the length scale of the turbulent field. 
When a is much larger than the Lorentz factor of a radiating electron 7, syn- 
chrotron radiation is realized, while a <^ 1 corresponds to the so-called jitter 
radiation regime. Because for 1 < a < 7 we cannot use either approximations, 
we should have recourse to the Lienard-Wiechert potential to evaluate the radia- 
tion spectrum, which is performed in this paper. We generate random magnetic 
fields assuming Kolmogorov turbulence, inject monoenergetic electrons, solve the 
equation of motion, and calculate the radiation spectrum. We perform numeri- 
cal calculations for several values of a with 7 = 10. We obtain various types of 
spectra ranging between jitter radiation and synchrotron radiation. For a ~ 7, 
the spectrum turns out to take a novel shape which has not been noticed up to 
now. It is like a synchrotron spectrum in the middle energy region, but in the 
low frequency region it is a broken power law and in the high frequency region 
an extra power law component appears beyond the synchrotron cutoff. We give 
a physical explanation of these features. 

Subject headings: magnetic fields — turbulence — radiation mechanisms: gen- 
eral — gamma-ray burst: general 

1. Introduction 

A major part of non-thermal emission from high energy astrophysical objects is almost 
always characterized by the radiation from relativistic electrons moving in magnetic fields. 



- 2 - 



Usually it is interpreted in terms of the synchrotron radiation. However, synchrotron ap- 
proximation is not always valid, in particular when the magnetic fields are highly turbulent. 
Electrons suffer from random accelerations and do not trace a helical trajectory. In general, 
the radiation spectrum is characterized by the strength parameter 

a = Ab 5", (1) 

where Ab is the typical scale of turbulent fields, \B\ is the mean value of the turbulent 
magnetic fields, e is the elementary charge, m is the mass of electron and c is the speed 
of light (Reville & Kirk 2010). When a ^ 7, where 7 is the Lorentz factor of radiating 
electron, the scale of turbulent fields is much larger than the Larmor radius = •jmc^ /e\B\, 
and electrons move in an approximately uniform field, so that the synchrotron approximation 
is valid. In contrast, when a <^ 1, Ab is much smaller than the scale rg/7 — X^/a which 
corresponds to the emission of the characteristic synchrotron frequency. In this regime, 
electrons move approximately straightly, and jitter approximation or the weak random field 
regime of diffusive synchrotron radiation (DSR) can be applied (Medvedev 2010, Fleishman 
and Urtiev 2010). For 1 < a < 7, no simple approximation of the radiation spectrum has 
been known. 

The standard model of Gamma Ray Bursts (GRB) is based on the synchrotron radia- 
tion from accelerated electrons at the internal shocks. The observational spectra of prompt 
emission of GRB can be well fitted by a broken power law spectrum which is called the 
Band function. Around a third of GRBs show a spectrum in the low energy side harder than 
the synchrotron theory predicts. To explain this, other radiation mechanisms are needed. 
Medvedev examined relativistic coUisionless shocks in relevance to internal shocks of GRB, 
and noticed the generation of small scale turbulent magnetic fields near the shock front 
(Medvedev h Loeb 1999). Then he calculated analytically radiation spectrum from elec- 
trons moving in small scale turbulent magnetic fields, to make a harder spectrum than the 
synchrotron radiation (Medvedev 2000). However, he assumed that the strength parame- 
ter a is much smaller than 1 and that turbulent field is of one-dimensional structure, which 
may be over simplified in general (Fleishman 2006). Medvedev also calculated 3-dimensional 
structure assuming that the turbulent field is highly anisotropic (Medvedev 2006). He con- 
clude that the harder spectrum is achieved in "head on" case, and that in "edge on" case, 
the spectrum is softer than synchrotron radiation. The spectral index depends on the angle 
6 between the particle velocity and shock normal with hard spectrum obtained when 9 < 10° 
(Medvedev 2009). 

Recently several particle-in-cell (PIC) simulations of relativistic coUisionless shocks have 
been performed to study the nature of turbulent magnetic fields which are generated near 
the shock front (e.g., Prederiksen et al. 2004; Kato 2005; Chang et al. 2008; HaugboUe 
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2010) . The characteristic scale of the magnetic fields is the order of skin depth as predicted 
by the analysis of Weibel instability. Then, the wavelength of turbulent magnetic field Ab is 
described by using a coefficient k as 

Ab - (2) 

l^pe^ int 

where is the plasma frequency, and Tint is the relative Lorentz factor between colliding 
shells. The energy conversion fraction into the magnetic fields 

is 10~^ — 0.1, where B'^/Stt is energy density of magnetic fields, and Fiut^mpC^ is the kinetic 
energy density of the shell. The Lorentz factor of electrons is similar to Tint, and that k is 
typically 10 from the result of PIC simulations, the strength parameter a can be estimated 
as 

a = ^Ke^^l\[^r^ 0(1 - 10^). (4) 
V rrie 

Thus, the assumption a <^ 1 on which jitter radiation and DSR weak random field regime 
are based is not necessarily valid when we consider the radiation from the internal shock 
region of GRB. 

Fleishman & Urtiev (2010) calculated the radiation spectrum for a > 1 using a statistical 
method, but their treatment of the "small scale component" is somewhat arbitrary. They 
introduced the critical wavelength Acrit and called components with A < Acrit the "small scale 
component" , where Acrit obeys the inequality 

Tg < Acrit < rg/7 (5) 

(Toptygin & Fleishman 1987). The inequality can be transformed to 1 ^ Acrite|-B|/mc^ <^ 7, 
so that the division through Acrit may be ambiguous when we calculate the radiation spectrum 
for 1 < a < 7. 

The synthetic spectra from PIC simulations were calculated recently (e.g., Hededal 
2005; Sironi & Spitkovsky 2009; Frederiksen 2010; Reville & Kirk 2010; Nishikawa et al 

2011) . Althogh their magnetic fields are realistic and self-consistent, it is inevitable that 
the fields are described by discrete cells in PIC simulation. Reville & Kirk (2010) developed 
an alternative method of calculation of radiation spectra that uses the concept of photon 
formation length, which costs much shorter time than the first principle method utilizing 
the Lienard-Wiechert potential. 

In this paper, we rather use the first principle method to obtain the spectrum as exact 
as possible. We adopt the field description method developed by Giacalone & Jokipii (1999) 
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and used by Reville & Kirk (2010). We assume isotropic turbulent magnetic fields which 
have broader power spectra k^^ax = 100 x /cmin and calculate the radiation spectra in the 
regime of 1 < a < 7. In §2 we describe calculation method and numerical results. In §3 we 
give a physical interpretation. 



2. Method of calculation 

Because we focus our attention on calculating radiation spectrum, we assume the static 
field with required properties of a, and neglect the back reaction of radiating electrons to the 
magnetic field. We solve the trajectory of electron accurately in each time step and calculate 
the radiation spectrum. 



2.1. Setting 

The isotropic turbulent field is generated by using the discrete Fourier transform de- 
scription as developed in Giacalone & Jokipii (1999). It is described as a superposition of A'^ 
Fourier modes, each with a random phase, direction and polarization 

N 

^i^) ^^^ri exp{i(fc„ • X + /3n)}in- (6) 
n=l 

Here, A„, Pn, fe„and^„ are the amplitude, phase, wave vector and polarization vector for the 
n-th mode, respectively. The polarization vector is determined by a single angle < < 27r 

in = COS i>ne'^ + i sin -0^6^, (7) 

where and are unit vectors, orthogonal to = kn/kn- The amphtude of each mode 



IS 

N -, -1 



EG" 



i-n=l 



(8) 



where the variance a represents the amplitude of the turbulent field. We use the following 
form for the power spectrum 

'^"-l + (fc„L,)«' 

where Lc is the correlation length of the field. Here, AA;„ is chosen such that there is 
an equal spacing in logarithmic A;— space, over the finite interval /cmin ^ k ^ /cmax, where 
^max = 100 X kjnin ^ud N — 100. whcrc /Cmin = Stt/Lc and a — 11/3. We have no reliable 



- 5 - 



constraint for value of a from GRB observation, so we adopt the Kolmogorov turbulence 
B'^(k) oc k~^^^ tentatively, where the power spectrum has a peak at fcmm- Then we define 
the strength parameter using a and fcmin as 

2'Kea 

a^^-—. (10) 

UlC rCnijjj 

We inject isotropically 32 monoenergetic electrons with 7 = 10 in the prescribed magnetic 
fields, and solve the equation of motion 

^rrieV = —e(3 x B (11) 

using 2nd order Runge-Kutta method. We pursue the orbit of electrons up to the time 
300 X Tg , where Tg is the gyro time Tg = 2Tr'ymc/ea. We calculate radiation spectrum using 
acceleration v = (3c. The energy dW emitted per unit solid angle dVt (around the direction 
n) and per unit frequency du to the direction n is computed as 



dW e 



2 



00 



dujdVL Attc^ 

where r{t') is the electron trajectory, t' is retarded time (Jackson 1999). 



2 



, ,n X \(n — f3) X 3] ,, n-r(t').. 

(1 — p ■ n)^ c 



2.2. Results 

First, we show the radiation spectrum for a = 3 in Figured] The frequency is normalized 
by the fundamental frequency Ug = ea/{'-ymc), and the magnitude is arbitrarily scaled. The 
jagged line is the calculated spectrum, while the straight line drawn in the low frequency 
region is a line fitted to a power law spectrum. The fitting is made in the range of 1 — 350a;g 
and the spectral index turns out to be 0.44. The straight line drawn in the high frequency 
region shows a spectrum of oc uj~^^^ expected for diffusive synchrotron radiation for reference 
(Toptygin & Fleishman 1987). The spectrum is well described by a broken power law, and 
the spectral index of the low energy side is harder than synchrotron theory predicts. The 
peak frequency of this spectrum is located at around lO^Ug. This frequency corresponds to 
approximately the typical frequency of synchrotron radiation Ugyn = 37^e(T/2mc ~ lO^Ug, 
for 7 = 10. 

Figure [2] shows the spectrum for a = 5. The spectral shape changes from that of a = 3 
in both sides of the peak. The spectrum of the low frequency side becomes a broken power 
law with a break around lOwg, above which the spectrum is fitted by a power law with an 
index of 0.33, as expected for synchrotron radiation, while below the break the index is 0.58. 
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The high frequency side above the peak indicates an excess above a power law spectrum 
dW/dudVL oc w"^/^. It looks like an exponential cutoff. The whole spectrum is described 
by a superposition of a synchrotron spectrum and a DSR broken power law spectrum. This 
spectral shape is totally novel and is different from the one by described Fleishman (2010). 
He reported that the spectrum is described by broken power law in the same range of a as 
this work (l<ci<7, ci~ 10*^ and 7 ~ 10^) (Fleishman 2010). 

To confirm our inference we calculate the case of a = 7, for three different values 
of a, i.e., a = 14/3, 11/3 and 8/3 and the results are shown in Figure O The curved 
black line is a theoretical curve of synchrotron radiation, and three straight black lines are 
dW/dcodil oc uj~'^^^,uj~^^^ and uj~^^^ expected for DSR theory for reference. The green line 
corresponding to a = 14/3 reveals a clear exponential cutoff, and reveals DSR component 
in only the highest frequency region. The power law index of this component in the highest 
frequency region coincides with the expected value —a + 2 = —8/3. The red and blue hues 
correspond to a = 11/3 and 8/3, respectively. They indicate a common feature to the green 
one. 



3. Interpretation 

We give a physical explanation of the spectra obtained in the previous section. At first, 
we interpret the broken power law spectrum for a = 3 by using DSR theory. Next, we 
consider physical interpretation of the complex shape of radiation spectrum for a = 5 and 
a = 7. Finally, we compare our work with previous studies. 

To begin with, we review the spectral feature of DSR based on the non-perturbative ap- 
proach for a < 1 (Fleishman 2006). The typical spectrum takes a following form: dW/dudQ oc 
w^/^ in the low frequency region, oc co^ in the middle frequency region, and oc a;"""''^ in the 
high frequency region. The low frequency part and the middle frequency part are separated 
at wlm ~ aUsyn- This spectral break corresponds to the break of the straight orbit approx- 
imation due to the effect of multiple scattering (Fleishman 2006). On the other hand, the 
middle frequency part and the high frequency part are separated at the typical frequency 
of jitter radiation uuu = t^jit which is estimated by using the method of virtual quanta as 
i^jit ~ 7^^minC ~ a~^Usyn (Medvedcv 2000, Rybicki & Lightman 1979). Then, for a ~ 1, 
wlm ~ i^MH is achieved and the middle region may vanish. The spectrum for a ~ 1 be- 
comes a broken power law with only one break, which is located at roughly the synchrotron 
frequency Ugyn- The power law index of low frequency side is ~ 0.5 (which is harder than 
synchrotron radiation), and that of high frequency side is —a + 2 = —5/3. Thus, the spectral 
feature for a = 3 can be explained by an extrapolation of DSR non-perturbative approach 
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for a < 1, if we consider that the middle frequency region is not conspicuous. Although 
our spectral index 0.44 slightly differs from 0.5 for DSR, this index is still harder than the 
synchrotron theory. The situation a ~ 1 can be achieved at the internal shock region of 
GRB, so that this may be responsible for harder spectral index than synchrotron observed 
for some GRBs. 

Next we interpret the spectral features for a = 5 and 7. The conceptual diagram of 
these spectra for 5 < a < 7 is depicted in Figure IH and a schematic picture of an electron 
trajectory is depicted in Figure El We explain the appearance of another break in the low 
frequency range seen at around lOcUg in Figure |2j On the scale smaller than Ab, the electron 
motion may be approximated by a helical orbit, while it is regarded as a randomly fluctuating 
trajectory when seen on scales larger than Ab. Therefore, for the former scale, we can apply 
the synchrotron approximation to the emitted radiation. The beaming cone corresponding 
to the frequency uj is given by 



(Jackson 1999). The deflection angle 9q of the electron orbit during a time Ab/c is estimated 
iohe 9q = a/'-) from the condition 



as seen in Figure [51 Thus, the synchrotron theory is applicable only for ^conc < ^O; so that 
the break frequency is determined by 6*0 = 6'cone5 and we obtain 



This break frequency is the same as obtained by Medvedev (Medvedev 2010). We understand 
that as a is larger, break frequency becomes lower, and when a is comparable to 7, Wbr 
coincides with the fundamental frequency eaj^mc. 

Next, we discuss on the high frequency radiation, which results from the electron tra- 
jectory on scales smaller than A^. The synchrotron theory applies between Ab/o = ''"g/7 
and Ab. However, we should notice that electron motion suffers from acceleration by mag- 
netic turbulence on scales smaller than X^ja. The trajectory down to the smallest scale of 
27r//Cinax IS jittering, which is attributed to higher wavenumber modes as seen in the zoom up 
of Figure |5l If the field in this regime is relatively weak, i.e., a is relatively large (Figure |3l 
green line: a = 14/3), the trajectory on the scale smaller than Ab/o does not much deviate 
from a helical orbit. In this case, radiation spectrum reveals an exponential cutoff, and a 
power law component appears only in the highest frequency region. On the contrary, if the 
smaller scale field is relatively strong, as in the case of a = 8/3 depicted in the blue line 
in Figure |3l the power law component becomes predominant in the high frequency region. 




(13) 
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and the synchrotron exponential cutoff is smeared out. The intersection frequency of curved 
black line and straight black lines at around lO^Wg in Figure [3] corresponds to wjit as seen 
in Figure HI Since the intersection frequency is determined by a, the frequency where the 
power law component appears over the synchrotron cutoff is dependent on a. The excess 
from the theoretical curve in the middle frequency region in Figure [3] may be explained by 
consideration of two effects. One is the contribution of hidden DSR component, and the 
other is a range of synchrotron peak frequency which is caused by a fluctuation of magnetic 
field intensity. 

Fleishman reported that the spectrum for 1 < a < 7 and 3 < a < 4 becomes a broken 
power law (Fleishman & Urtiev 2010). Medvedev asserted that the high frequency region 
of the spectrum reveals an exponential cutoff for 1 < a < 7 (Medvedev 2010). Our result 
indicates that an exponential cutoff plus an extra power law component appears, which 
is different from Fleishman's remark and from Medvedev's remark on the high frequency 
region. On the other hand, similar spectra to ours have been reported in Fleishman (2005) 
and Reville & Kirk (2010) when a uniform field is added to turbulent field. Because the high 
energy power law component arises from a turbulent spectrum over the wavenumber space, 
this component does not exist when the small scale field is excited only in a narrow range 
of wavenumber space. Since the energy cascade of turbulent magnetic fields should exist 
at least to some degree, we regard that the higher wavenumber modes naturally exist. It 
depends on the set of parameters of crandfcmax whether this high energy power law component 
can be seen or not. If 27recr/mc^A;inax > 1? this component will not be seen. If the magnetic 
turbulence is excited by Weibel instability at the relativistic shocks, it is not possible for 
^max to be much larger than fcmin because the wavelength of injection (Ab = I-k jk^y^ is only 
a few ten times the skin depth at most. Therefore, the component will not be seen for a 3> 1 
while for a ~ 0(1), this power law component will be seen. 

As for the frequency region lower than the break frequency Wbr = '^"^i^syn; Medvedev 
remarked that the spectrum is similar to small angle jitter radiation (Medvedev 2010). How- 
ever, it remains to be open if it is so for 1 < a < 7, because the assumption that the straight 
orbit approximation of radiating particle is broken. To predict the exact radiation spectrum 
of the frequency region lower than the break frequency, it is necessary to pursue the particle 
orbit to follow the long term diffusion which is a formidable task. 

4. Summary 

We calculate the radiation spectrum from relativistic electrons moving in the small 
scale turbulent magnetic fields by using the first principle calculation utilizing the Lienard- 
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Wiechert potential. We concentrate our calculation on a range of the strength parameter of 
1 < a < 7. We confirm that the spectrum for a ~ 3 is a broken power law with an index 
of low energy side ~ 0.5, and that some GRBs with low energy spectral index harder than 
synchrotron theory predicts may be explained. Furthermore, we find that the spectrum for 
a ~ 7 takes a novel shape described by a superposition of a broken power law spectrum and a 
synchrotron one. Especially, an extra power law component appears beyond the synchrotron 
cutoff in the high frequency region reflecting magnetic field fluctuation spectrum. This is 
in contrast with previous works (Fleishman & Urtiev 2010, Medvedev 2010). Our spectra 
for a = 5 and a = 7 are different from both of them. We have given a physical reason 
for this spectral feature. This novel spectral shape may be seen in various other scenes in 
astrophysics. For example, the spectrum of 3C273 jet at the knot region may be due to this 
feature (Uchiyama et al. 2006). 

We thank the referee for helpful comments. We are grateful to T. Okada, S. Tanaka, M. 
Yamaguchi for discussion and suggestions. This work is partially supported by KAKENHI 
20540231 (F.T.). 
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Fig. 1. — Radiation spectrum for a = 3, a = 11/3 and 7 = 10. The straight hne in the low 
frequency region shows a power law spectrum with an index 0.44. The straight line in the 
high frequency region is dW/dco cx cu~^^^ for reference. Power law index of the low frequency 
spectrum is harder than the synchrotron theory predicts. 
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Fig. 2. — Radiation spectrum for a = 5, a — 11/3 and 7 = 10. Two straight lines drawn in 
low frequency region are fitted one to a power law spectrum in the range of 0.5 — lOug and 
10 — lO'^cjg, respectively. The power law index is 0.58 and 0.33. The latter corresponds to the 
synchrotron radiation. The straight line in the high frequency region is dW/dcu oc uj~^^^, and 
we see a broad hump in the peak region, which is identified with the synchrotron spectrum 
with an exponential cutoff. 
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CO/COg 

Fig. 3. — Radiation spectra for a = 7 and 7 = 10, with a = 14/3, 11/3 and 8/3. Green, red 
and blue lines are calculated spectra a = 14/3, 11/3 and 8/3, respectively. Straight black 
lines are dW/du oc o;"""'"^. Curved black line is a theoretical curve of synchrotron radiation. 
We see effects of different spectrum of the turbulent magnetic fields in high energy region. 
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Fig. 4. — Conceptual diagram of the radiation spectrum for 5 < a < 7. 
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Fig. 5. — Cartoon showing the electron trajectory for 1 < a < 7. The radius of the guiding 
circle is the Larmor radius = 'jX^/a. Low frequency photons are emitted from the motion 
on relatively scales larger than Ab, while middle frequency ones from that on the intermediate 
scale between Xb/o, and Ab are basically synchrotron radiation. The spectral break at a'^j^ 
in Figure m corresponds to the break of synchrotron approximation at the scale of Ab- The 
scale Ab/a = f^g/l corresponds to the synchrotron peak frequency. On the smallest scale 
down to 27r/A;max; the trajectory is approximately straight, and jittering is responsible for 
the power law component in the highest frequency region in Figure |H 



